Analogy of RKKY oscillations to the heat exchange in cold atoms 
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An oscillatory term is found in both the energy expectation and dynamics of a wave-packet in a 
time-varying harmonic trap and infinite potential well. They are proved to oscillate in coherence 
with the time lapse within each period depending on both the cutoff in transition energies and the 
specific route via which the potential is being varied. This oscillatory term is general to arbitrary 
potential forms since it derives from the interference between crossed transition trajectories. Close 
analogy is made to the Ruderman-Kittel-Kasuya-Yosida interaction for giant-magnetoresistance 
trilayers, where many-body quantum interference among scattering states renders the oscillation 
as a function of spacer width. This connection reveals the generality of quantum friction due to 
parasitic oscillations. 
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I. INTRODUCTION 

The adiabatic and sudden processes are two opposite 
limits in the time-dependent quantum mechanics. The 
adiabatic tuning is involved in experiments, such as the 
advanced cooling technique[l[, quantum computing0,Q, 
and the state broadening in optical lattices [J]. The com- 
bination of adiabatic and sudden processes is also used 
to minimize the transition time between different ther- 
mal equilibrium states Q. Although analytic forms con- 
veniently exist for these two limits, cares need to be taken 
when identifying them to real systems @. How to de- 
scribe the physics for a time-varying potential when nei- 
ther limit is applicable remains a challenging task. In- 
teresting phenomena have been proposed to exist in this 
regime, for instance, the parasitic oscillations that give 
rise to the quantum friction when the Hamiltonians at 
different time do not commute, as is the case when the 
potential profile changes at finite rates 0, 

Unlike the friction that lead to dissipations, the quan- 
tum frictional process is reversible and energies can 
be stored in or taken out of the oscillations Non- 
adiabatic processes allow the transition between differ- 
ent instantaneous states, which leads to the parasitic os- 
cillations. Ostensibly it seems difficult to deduce any 
general property in these systems without the knowledge 
of the specific paths of transition. However, since the 
parasitic oscillations collective behavior that re- 

sults from the contribution of all trajectories, this task 
has been shown to be possible. One of the most famous 
examples in many-body systems is the carrier-mediated 
interaction between magnetic impurities or thin films, 
namely, the Ruderman-Kittel-Kasuya-Yosida (RKKY) 
coupling [9]. This coupling is dominated by the quan- 
tum interference between scattering paths [10( of the me- 
diating electrons and the resulting strength oscillates as 
a function of separating distance with a period of half 
Fermi wavelength. Similar to the RKKY coupling, the 
dynamic parasitic oscillations mediate the coupling be- 
tween boundaries. While the energy input from the mov- 
ing boundaries plays the role of coupling strength, the 
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FIG. 1: (color online) Interference diagrams for (a) the RKKY 
oscillation and (b) oscillatory term in the energy expectation 
in a contracting potential well. The curve lines denote scatter- 
ing paths with different rounds between the interfaces. The 
solid lines represent energy levels and the dotted lines are 
transition paths. 



oscillation is now in time instead of space. If we take 
as an example the giant magnetoresistance (GMR)fll| 
system in which the RKKY mechanism has successfully 
applied, the close resemblance of interference diagrams 
for RKKY and the time-varying potential trap is exem- 
plified in FigHJ 

In Section II we first examine the harmonic trap be- 
cause the time-dependence of its basis lies solely in the 
frequency oo(t). The expectation values of energy and 
position-momentum correlation that represent the para- 
sitic oscillations are calculated analytically for an arbi- 
trary wave-packet. We extend our study in Section III to 
an infinite well to demonstrate that these oscillations and 
the time lapse within each period they share is calculable 
are general. The interesting analogy between RKKY os- 
cillation and the energy input from the time- varying trap 
is arranged in the final Section IV, together with some 
discussions and conclusions. By use of formulae in Sec- 
tion II, we calculate the energy expectation of harmonic 
oscillators varying with a constant adiabatic parameter 
in Appendix A, and reproduced routes that have been 
reportedpjj to be frictionless. Appendix B provides cal- 
culations to explain why the oscillation ceases to exist in 
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the adiabatic and sudden limits. 



II. HARMONIC TRAP 

A. Energy and Position-Momentum Correlation 

We consider an one-dimensional time-dependent har- 
monic trap, H(t) = p 2 /2m + mui 2 (t)x 2 /2, with an 
initial angular frequency w(0) = loq. An orthonor- 
mal basis can be constructed out of the time-dependent 
cigenfunctions (l3j 
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where b{2h/niLOo) 1 / 2 = bLo is the time-dependent length 
scale with the scaling function, b(t), satisfying 



b + UJ 2 {t)b = UJ^/b\ 



(2) 



and the initial conditions 6(0) = 1. 

To calculate the energy expectation, we first use H(t) 
to operate on the basis function 

H(t)\^ n (t)) = e~ i{n+1/2) ti dt '("o/b 2 ) e i(7n/2h)bx 2 /b 

(n + \P£ + j(b 2 bb)(* b ) 2 + ±{x,p}] |0 n (6)), 
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where \4> n (b)} is the normalized instantaneous eigen- 
function for {p 2 /2m + muj 2 x 2 /2b ) with eigenvalue (n + 
l/2){hw Q /b 2 ). 

Given any wave-packet ip(0,x) = ipo(x), its evolution 
follows ip(t,x) = J2n c n^n(t,x) with the projection 

cn = (*n(* = 0)1-00) = (<t>n{b = l)\e- l(m / 2h ^ x2 ip ). 

(4) 

By use of Eq.®, the energy expectation can be obtained 

as 
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Since our choice of time-dependent basis is not the eigen- 
functions of energy operator, the energy expectation 
in Eq.([5|) contains both diagonal and non-diagonal el- 
ements. The non-diagonal elements are non-zero only 
when n and / share the same parity; namely, n — I is an 
even number. 



The two summations in Eq.([5j) can be removed by 
appealing to the identity operator after we change the 
length scale in the second inner product from b to 1 and 
reexpress the (n — l) in the exponent by the Hamiltonian 
operator. Equation ((5|) now becomes 

(H(t)) ^\e mh ^' l2 / 2h e 1 ^ /o dt '/ b2 

[(fs*.)^,^] 

with Ho = p 2 /2m + muJ^x 1 /2. Being independent of b, 
the first term in the square bracket can be identified as 
the energy in the adiabatic limit: 
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where Eq.([2]) has been used. In the mean time, the 
mb 2 x 2 /2 term describes the extra kinetic energy due to 
the pull or push by the moving potential. The third term 
can be easily shown to be proportional to [H(0), H(t)], 
which gives rise to the quantum friction 0, H[ when 
the Hamiltonians at different times do not commute. 
The reason is that the positive/negative sign of ({x,p}) 
reflects the expansion/contraction motion of the wave- 
packet, while b/b keeps track of similar motions of the 
trap. If the wave-packet moves along with the trap, it 
receives a positive work of b({x,p}) /2b. On the other 
hand, against intuitions, the packet can sometimes resist 
and behave oppositely, which causes the work to turn 
negative and diminishes its energy. If the initial wave- 
function consists of only one ipo Eq.([T|) as, Eq.© can be 
solved straightforwardly as 
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The friction term expectedly becomes zero because the 
eigenfunction always moves coherently with the trap. 

For a general ipo, Eq.© can be simplified by using (1) 
ladder operators in Hq, ao and aj, and (2) the commu- 
tation relations 
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By use of both relations, Eq.© can be simplified to 
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where the new ladder operators, A and A + , are obtained 
from the old set via the Bogoliubov transformation, 
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Another interesting property pertinent to the dynamics 
is the expectation value for the position-momentum cor- 
relation operator C = {x,p}, 
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The derivations so far depend on our selection of a time- 
varying trap that avoids generating a complex phase fac- 



tor, ujq = \J cj 2 (0) + 6(0). One example is the trap with 
a constant adiabatic parameter /i = lu/uj 2 . The energy 
expectation was found to switch from being sinusoidal to 
hyperbolic as \fj,\ > 2 by Rezek et al.\T^. Furthermore, 
they discovered the oscillation to be intimately linked to 
the existence of frictionless routes when < 2. Detailed 
discussion are arranged in Appendix A where identical 
results are rederived by our approach, i.e., Eq.®. From 
now on, we shall confine ourselves to initial conditions 
6(0) = that guarantees u>o = ui(0) to be real. 

It may seem surprising for a system as simple as the 
harmonic potential within the single-particle picture that 
the energy expectation in Eq.® should contain an oscil- 
latory term which, within one period, requires 



dt'uj /b' 2 (t) = 2tt. 
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where T is in general time-dependent and thus is not a 
proper period. To understand its origin, let us go back 
and examine the mathematic structure in Eq.([5]). We 
have already shown that only the even-(n — I) terms need 
to be considered. Therefore, the T shared by Eq.© and 
Ea. ([TT1) is the least common multiple of these different 
phases as defined by Eq. (fT7|) . 



B. Coupling between quantum friction and 
quantum breathing 

The parasitic oscillations give rise to the breathing 
effect [1J] which term was originally proposed to describe 
the large-amplitude fluctuations of wave-packet caused 
by the excess energy [l5j when a harmonic trap was ex- 
panded suddenly. In this case, the energy expectation 
of the final state is a constant. Only the position- 
momentum correlation (C) contains the oscillatory term 
with frequency, 1/T, determined solely by the frequency 
u> characterizing the final trap. In the following, we shall 
allow the trap to contact at a general rate. For a medium 
rate (66/wo) ~ 1, the quantum friction will enter and ren- 
der the quantum breathing with the oscillation period 



sensitively affected by the specific route via which the 
final stage is reached. 

We start with a linearly contracting harmonic trap, 
where 6 is constant, to pave ways for our later extension 
to other potential profiles in Section III. Set the initial 
wavefunction to lie in the ground state \4>o(b = 1)) and 
insert it to Eq.® give the energy expectation at any later 
time 
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For a mild contraction, (66/wo) ~ 1, the oscillation terms 
are as important as the first term in Ea. (fl"3|) . Besides 
demonstrating the oscillatory feature, FigJ5]also shows a 
higher energy expectation than the adiabatic result which 
is an artifact of our choice to start from the ground state. 
Had we chosen a mixed state such as the coherent state 
cxp(ipoi/^) exp(— ixop/h) |</>o(l)), where xo/po denotes 
the initial position/momentum of the state, their relative 
size could be reversed. For the coherent state that moves 
toward the contracting boundary, the energy modulation 
is smaller than the adiabatic result because of the reduc- 
tion from the quantum friction. 

Another useful signature of the wave-packet dynamics 
is the position-momentum correlation (C) in Eq. (|lll) . 
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where C(0) = and a positive/negative C(t) signifies 
an expanding/ contracting wave-packet. The numerical 
result in Fig(3] reveals an unexpected motion; namely, 
the wave-packet can sometimes expand against the con- 
tracting trap. In FiglU the variance var(x)= (a: 2 ) — (x) 2 
is calculated and plotted in the solid line. Consistent 
with the conclusion of Figj3l the uncertainty in position 
also increases in the gray regions, in clear contrast to the 
monotonous adiabatic result in the dashed line. 

Since our choice of initial wavefunction (f>o(b = 1) has 
nonzero projections on more than one direction of the or- 
thonormal basis & n (t, x) in Eq.([T]), the quantum breath- 
ing or the parasitic oscillation of its density profile is ex- 
pected. Similar to its classical counterpart, the quantum- 
friction term extracts work from and thus slows down 
the wave-packet when its motion is against that of the 
potential profile. This resistance will modulate the en- 
ergy expectation as well as its dynamics and result in 
oscillations coherently but not necessarily synchronized. 
This is demonstrated by the contrast of the solid and 
dashed lines in Fig|3J which monitor the dynamics and 
the amount of energy input, respectively. 
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FIG. 2: (color online) The energy expectation in Eq. (|13p 
is plotted as a function of b for a linear contraction with 
b = — a>o/4 and b — (thick solid line), and the adiabatic 
process with E(b) = £t)/6 2 (thin solid line). In compari- 
son, the coherent state with {xQ,po}={^/h/mu)o, — y/hmujo} 
that moves away from the contracting boundary is plotted 
in dashed line and the state with the opposite motion for 
{xo,po}={y /, h/mujo, Vftmojo} is plotted in dotted line. From 
now on, the energies in all figures are expressed in unit of the 
initial energy. 




FIG. 3: (color online) The solid line denotes the expectation 
value of {x,p} in Eq. (|14[) with the same parameters as its 
counterpart in Fig[2] For comparison, the difference in en- 
ergy expectation between the two lines in Figgis plotted in 
the dashed line. The gray areas highlight the regions when 
the wave-packet expands against the contraction of the trap. 
Although the energy input remains positive when this unex- 
pected behavior occurs, its efficiency is greatly diminished due 
to the quantum friction. 



However, the harmonic oscillator is special because of 
its unique internal frequency and the classical analogies 
it enjoys, such as the dynamics of the coherent state. In 
order to demonstrate (1) the authenticity of this oscil- 
latory feature, (2) the dynamics and energy are coupled 
together with a common time lapse that is determinable, 
and (3) our approach of constructing a basis out of the 
original cigenfunctions is general to all potential traps, 
we shall examine the infinite potential well in the next 
section. 



FIG. 4: (color online) The variance, var(a)= {x 2 } — (a;) 2 , is 
plotted as a function of b for b — — ojo/4 (solid line), — uo/2000 
(dashed line and adiabatic limit), and — 25ojo (dotted line and 
sudden limit). The length scale Lo is defined as (2K/mwo) • 
The gray area are the same regions of interest in FigO 



III. INFINITE WELL 

For any one-dimensional linearly-varying trap, the 
hamiltonian can be described by H(x/(bL )) with the di- 
mensionless factor b{t) obeying 6 ^constant and 6(0) = 1. 
One convenient orthonormal basis can be chosen as 



,(t,x) =b-ie- i ^dt'E n /H e imb^/2bfl ( f ) JX A (1 5 ) 

\bLnJ 



where n are positive integers and 4>n( b ^ Lo ) and E n (t) are 
the instantaneous eigenstates and eigenencrgies. The ex- 
ponents, (— f Q dt'E n (t')/h) and (mbx 2 /26ft), provide the 
dynamic phases. By using this basis, the energy expec- 
tation for any initial state ?po can be found to share the 
same structure as Eq.© with 6 = 0. Therefore, the ar- 
guments following Eq.Q can also be applied here, which 
allow us to conclude that the coupling between quantum 
friction and quantum breathing is general. 

For illustration, let us pick an initial wave-packet in 
Boltzmann distribution of the instantaneous eigenstates 
at t = 0: 



Mx)=N- 1 / 2 J2 e ~ (3Eni0)( t>n(x/L ), 



(16) 



n=l 



where N = J^exp(— /3E n (0)). Both the oscillation fea- 
ture and the coupling between dynamic and energy are 
displayed clearly in the numerical result in Fig JS] for the 
operating temperature ~ IK accessible to cold-atom ex- 
periments. Naively, one may not expect to see oscillations 
since there is no intrinsic frequency scale for an infinite 
well. It turns out that the phase (n — I) f Q dt'(ujo/b 2 ) in 
Eq.([5]) can be reexpressed as J ' dt'(E n — Ei)/H. And this 

gives rise to a cutoff, f Q dt'(Ei — E2)/fr, which defines a 
time lapse for each period 
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FIG. 5: (color online)The solid line denotes the expectation 
value of {x,p} with the b = #i(0)/2/i, (3 = l/2E 1 (0). As in 
FigEl the dashed line represents the energy expectation after 
subtracting the adiabatic contribution. 

This was checked to be consistent with the numerical 
results in Fig|5] It comes as no surprise that we also 
find the above mentioned features persist in the limit 
of f3 — > oo; namely, when the initial wave function is a 
pure state. In the mean time, we checked that the same 
oscillatory feature in the energy expectation was retained 
after the thermal averaging within the density matrix 
method. 

For fixed initial and final boundary positions, we now 
examine the efficiency of different contracting rates for 
the energy input. Figure [5] shows the extra energy input 
compared to the adiabatic contribution oscillates with 
the contracting time. In contrast to the similar definition 
in Eq. (|17p , the time lapse T is now a constant for a fixed 
b and thus is a well-defined period. Since our derivations 
are applicable to all trapping systems by the choice of 
basis in Eq.(|15p. we argue that the temporal oscillation 
in energy expectation with a fixed period is a common 
feature for all linearly-varying traps. 

IV. DISCUSSIONS AND CONCLUSIONS 

The moving boundaries are coupled through the atom 
trapped inside, and the energy input can be used to define 
the strength of this coupling. The oscillation with con- 
tracting time is reminiscent of the dependence on spatial 
distance between magnetic impurities in the Ruderman- 
Kittel-Kasuya-Yosida (RKKY) interaction mediated by 
the conducting electrons. The coupling between mag- 
netic side-layers in trilayer systems, such as the giant- 
magnetoresistence samples, has been shown [l6j to be a 
simple extension of the RKKY formula: 

v / n— 1 

(18) 

where A is the interface area, Ep is Fermi energy, Tr/ti, 
represent the refraction coefficients in right/left inter- 



FIG. 6: (color online)The extra energy input in addition to 
the adiabatic result is shown for different contracting times 
in the unit to = 2h/Ei(0). The final boundary is chosen 
arbitrarily to be at b = 0.3 for the solid line and 6 = 0.2 for the 
dashed line. Energies for the latter case have been suppressed 
by 10 times in order to fit in the plot. A short t signifies a 
fast contraction. This oscillatory feature is reminiscent of 
that of RKKY coupling when plotted as a function of spatial 
distance. 



Oscillation 


Energy Input 


RKKY coupling 


Definition 
Ensemble 
Contributors 
Parameter space 
Period 


E(t) - E(t -> oo) 

Canonical 
Transition paths 
Contracting time 
hb/E c (0) 


E(D) - E(D oo) 
Grand Canonical 
Scattering paths 
Spatial distance 
n/ftF 



TABLE I: Comparison of quantum oscillations in the energy 
input for a time- varying trap and in the RKKY coupling. 



faces, and D denotes the spacer width. The perpendicu- 
lar wavevector k± is a function of energy E and parallel 
wavevector . The above equation indicates that the 
coupling energy is contributed by all possible scattering 
paths below the Fermi sea. Although Eq. (fT5|) can not 
be solved analytically, numerical results show that it is 
an oscillatory function in D with period ir/kp which is 
half of the Fermi wavelength. This period is expected 
because (1) the integrand consists of sinusoidal functions 
and (2) the Fermi energy defines a natural cutoff in the 
period[10|. In comparison, the oscillatory energy input 
in the time-varying traps is the result of quantum in- 
terference between different transition paths, and it also 
contains a well-defined period as T = hb/E c (0), as shown 
in Eq. (fTT|) . Here E c (0) plays the role of energy cutoff as 
the Fermi energy in the RKKY scenario. Table U details 
the comparisons between these two oscillations. 

In conclusion, we find that the energy expectation in 
time-varying traps is coupled to the dynamic fluctua- 
tions through quantum frictions. Both the energy in- 
put and position-momentum correlation exhibit oscilla- 
tory features which we ascribe to the quantum interfer- 
ence among different transition paths. The time lapse in 
each period can be determined by the cutoff in transition 
energies. By fixing the initial and final trap profiles, we 
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find that the energy input is also an oscillatory function of 
the contracting time with an universal period hb/E c (0) 
for all linearly-varying trap. Mathematically we show 
that the formalism of this problem is analogous to that 
of indirect magnetic coupling in metals. Therefore, the 
oscillation in the time domain of the time- varying trap 
shares the same physical origin as the RKKY oscillation 
in the spatial domain. 



Appendix A: Harmonic Oscillator with a constant 
Adiabatic Parameter 



For a harmonic oscillator varying with a constant adia- 
batic parameter // = Cj/ut, the timc-dcpcndcnt frequency 
is 



;(«) = 



w(0) 



1 - (j, ■ u(0)t' 



(Al) 



The phase factor can be solved as ujq = uj(0)y/l — n 2 /4: 
by Eq.©. It is obvious that the phase factor is only 
real for < 2. The energy expectation is calculated via 
Eq.© as 



{ H(t)) =^{[i + ^ |A+A|^ ) 



-2/ 



^ t >W2 + * 7 )«>o|AA|V>o)]}, 

(A2) 

where 7 = ^1- (i 2 /4 and 6(t) = ln(^)//i. The ra- 
tio between the oscillatory and non-oscillatory term is 
time- independent, which reveals that the energy expec- 
tation reduces to the adiabatic result when the time lapse 
passes each period of the oscillatory term. For example, 
if we can select the ground state for the frequency w(0) 
as the initial wave-packet, the energy expectation can be 
derived from Eq. (|A2l) 



E(t) 



2 7 2 



;(2 7 0(t)) 



(A3) 



The condition 2j8(t n ) = 2mr defines the time lapses at 
which the energy returns to the adiabatic results ftiu(t)/2 
or, in other words, becomes frictionless. These observa- 
tions are consistent with those reported in Ref. 12]. Our 
Eq. ([A3p is equivalent to their Eq.(8) after replacing the 
2 7 by ifl. In contrast to our approach that starts from 
the initial wavefuntion, Salamon and Rezek[f| [l2[ esti- 
mated this energy by connecting the expectaion values 
of Hamiltonian H, Lagrangian L and correlation C for a 
time-dependent harmonic oscillator. It will be heuristic 
to clarify the transition at fi = 2 by comparing our for- 
mulae, Eqs.© and (fTTj) , and Salamon and Rczek's in a 
future research. 



Appendix B: Adiabatic and Sudden limit 

To get a more comprehensive understanding of the os- 
cillatory feature in time- varying traps, let us examine 
Eq.© in the adiabatic and sudden limits. The compe- 
tition between the time scales of the wave-packet and 
the trap variation is best exemplified in the special case 
of a linear expansion or contraction, b = v/Lq and 
b = 0. The adiabatic condition can be imposed by ei- 
ther (1) \\/2lo / {Suj 2 )\ -C 1, which states that the en- 
ergy input from the time-varying trap is much smaller 
than the intrinsic energy of the wave-packet (l7j or the 
more direct (2) |w/w p | -C 1 which restricts the potential 
to varying much slower than the averaged atomic speed, 
v p = \p2hj (mbLo), in a trap with an intrinsic angular fre- 
quency, uj(t) = LUo/b 2 . Both approaches guarantee that 
(bb/uJo) -C 1 and Eq.© can be simplified to give the 
energy in the adiabatic approximation as 

<*e»-£(" + i)-^?», (BD 

where N = (a + a). The oscillation term can be neglected 
in this limit. 

The sudden approximation describes the other limit of 
(bb/u)o) 3> 1 when the atomic motion can not catch up 
with the trap variation. The Taylor expansion works for 
the exponential term in Eq.© at this limit 



-2if*dt'u /b 2 
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Keep the expansion up to the fourth order and insert to 
Eq.© while retaining the lowest order term in LJo/(bb) 
give the energy expectation as 



(H) 



M^+i) + ^(^-l)((a + a + ) 2 ) 



(B3) 



with ui(t) = u>o/b 2 . This is the same result as is ob- 
tained by the time-independent perturbation theory. It 
is worth mentioning that, although this conclusion is in- 
dependent of the sign of b for a harmonic potential, the 
sudden approximation will break down terribly when ap- 
plied to a contracting infinite well because of the shrink- 
age of Hilbert space. Strictly speaking, the oscillatory 
term that we highlighted in Eq.© does not show up in 
the adiabatic and sudden limits for different reasons. In 
the former, it is suppressed due to the smallness of the co- 
efficient in Eq. © , while the Taylor expansion of Eq. (|B2p 
guarantees its minuteness in the latter. Finally, note that 
the energy expectation does not depend on b explicitly 
in these two limits. 
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